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ABSTRACT 


A  wave  height  transfonnation  model  is  developed  using  surface  roller  theory.  Roller  energy  production 
is  included  in  the  energy  flux  balance  equation  to  predict  rms  wave  height  for  randomly  varying,  irregular 
waves  over  arbitrary  bathymetry.  The  dissipation  function  is  defined  using  wave  roller  theory,  where  the 
area  of  the  roller  is  defined  from  a  simple  bore  analogy.  The  Rayleigh  distribution  is  used  to  statistically 
describe  wave  heights  as  waves  shoal,  break,  and  dissip^e.  Model  predictions  are  compared  with  data 
acquired  on  both  barred  and  near  planar  beaches.  The  sur&ce  roller  wave  height  transformation  model 
predicts  rms  wave  heights  with  an  average  rms  error  of  6.5%  for  a  barred  beach  over  three  days,  3.0% 
for  two  planar  beaches  over  four  days,  and  within  4.5%  average  error  for  all  locations  over  seven  days. 
The  model  has  two  free  parameters,  o  representing  the  type  of  breaker,  and  y  a  measure  of  breaking  wave 
saturation,  also  a  function  of  beach  slope.  Optimal  values  of  both  parameters  are  chosen  by  model  fitting. 
The  model  is  sensitive  to  y,  but  not  o.  The  sur^e  roller  model  improves  the  bore  dissipation  model 
[Thornton  and  Guza,  1983]  by  decreasing  the  average  rms  error  by  40%  while  decreasing  model 
sensitivity  to  input  parameters. 
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I.  INTRODUCTION 


Wave  heights  gradually  increase  as  waves  shoal  from  deeper  to  shallower  water.  As  waves  propagate 
into  very  shallow  water,  they  become  unstable  and  break,  generating  turbulence  at  the  sur^e  boundary 
layer.  The  turbulence  due  to  wave  breaking  is  the  primary  dissipative  mechanism  causing  wave  energy 
decay  in  the  surf  zone.  Dissipation  rate  and  turbulent  penetration  depth  are  dictated  by  breaker  type, 
commonly  described  as  spilling,  plunging,  or  surging.  While  plunging  wave  turbulence  can  penetrate  to 
the  bottom,  spilling  breaker  turbulence  is  primarily  confined  to  the  sur&ce  layer  between  the  crest  and 
trough.  Two  important  parameters  which  distinguish  breaker  type  are  wave  steepness  and  beach  slope. 

Cacina,  [1989]  concluded  that  wave  breaking  is  a  function  of  both  deep  water  wave  steepness  (Hq/Lo) 
and  beach  slope  (tan3),  and  is  correlated  with  the  deep  water  surf  similarity  parameter  described  by 
Battjes  [1974]  as 


C 


tanp 


H. 


miO 


N  h 


(1) 


where  is  the  deep  water  wave  height,  and  L„  is  the  deep  water  wave  length  (gTV27i). 

Several  methods  have  previously  been  employed  to  describe  wave  height  distributions  in  the  surf  zone. 
The  earliest  models  describe  shoaling  waves  as  monochromatic  and  entirely  dependent  on  local  water 
dq)th 
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where  y  is  an  adjustable  coefficient,  and  d  is  local  water  de(^.  More  recent  wave  height  tiamfoimation 
models  use  the  concept  of  energy  flux  balance  and  a  physically  based  dissipation  function.  In  the 
stochastic  models  of  Battjes  and  Janssen  [1978]  and  Thornton  and  Guza  [1983]  (TG83),  wave  heights  are 
described  statistically.  Energy  dissipation  due  to  shallow  water  wave-breaking  is  modeled  with  simple 
periodic  bores. 

The  objective  of  this  paper  is  threefold:  to  improve  the  previously  developed  bore  dissipation  model 
by  Thornton  and  Guza,  [1983]  by  incorporating  sur^e  roller  theory;  to  determine  the  influence  of  wave 
rollers  on  wave  height  transformation  in  the  surf  zone;  and  to  reduce  sensitivity  to  quasi-arbitiary  model 
parameters.  Waves  are  modeled  as  random,  having  irregular  amplitude,  and  approaching  from  angles  over 
arbitrary  bottom  profiles.  Roller  model  generated  wave  heights  are  compared  with  observed  wave  heights 
acquired  on  a  naturally  barred  beach  during  DELILAH  (Duck,  N.C.,  1990),  and  on  near-planar  beaches 
from  NSTS  (Leadbetter  Beach,  Santa  Barbara,  1978  and  Toney  Pines  Beach,  San  Diego,  California, 
1980). 
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II.  SURFACE  ROLLER  MODEL 


As  the  wave  shoals,  its  height  increases  and  the  wave  eventually  breaks  generating  turbulent  kinetic 
energy  (tke)  in  a  volume  of  water  which  becomes  detached  (separated)  from  the  wave  form  and  pitched 
forward  down  the  wave  fiace,  defined  as  the  wave  roller.  The  surface  roller  is  believed  to  play  a  vital  role 
in  wave  height  transformation  within  the  surf  zone  [Fredsoe  and  Deigaard.  1992].  Svendsen  and  Madsen 
[1984]  first  modeled  breaking  waves  as  rollers,  in  which  they  envisioned  a  mass  of  water  pushed  by  the 
wave  front  with  horizontal  velocity  equal  to  the  wave  celerity  C  (Fig  1).  The  shear  layer  between  the 
water  in  the  roller  and  wave  motion  results  in  the  dissipation  of  roller  tke  locally  within  the  wave/roller 
system.  Because  the  sur&ce  roller  propagates  with  the  wave,  roller  turbulence  can  be  the  main  source  of 
turbulence  in  the  inner  surf  zone  region  [Basco,  1983]. 

The  energy  flux  balance  equation  is  used  to  determine  wave  energy  gradients  across  the  surf  zone 


4<ECjsom)  =  D  (3) 

at  * 


where  E  is  energy,  is  the  wave  group  velocity,  a  is  the  mean  incident  wave  angle  calculated  using 
Snell's  Law,  and  D  is  dissipation.  As  in  TG83,  the  model  assumes  stationary  wave  conditions,  straight  and 
parallel  bottom  contours,  and  random  waves  which  are  narrow-banded  in  both  frequency  and  direction 
(i.e.,  all  waves  are  from  the  same  direction  and  of  a  single  frequency).  Energy  flux  is  contained  in  two 
terms,  rqiresenting  contributions  from  the  wave,  and  the  roller,  E„ 


E  ^  E  +  £ 


(4) 
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The  dissipation  term  in  Equation  ( 1)  is  derived  from  the  wave  roll«^.  Thus  Equation  (3)  becomes 


4:(EjC/x»a)  *  =  e,  (5) 

dr  *  dx 

where  e,  represents  roller  dissipation. 

Thornton  and  Guza  [1983]  show  that  wave  height  probability  density  functions  both  inside  and  outside 
the  surf  zone  are  well  described  by  the  Rayleigh  distribution,  p(H),  where 


piH)  =  -^«p[-(-^)^ 
HL.  "rm 


(6) 


The  ensemble  averaged  wave  energy,  <E^>,  is  calculated  by  integrating  through  the  Rayleigh 
distribution  using  linear  wave  theory 


(7) 


where  p  is  seawater  density,  g  is  gravitational  acceleration,  H  is  wave  height,  and  is  the  root-mean- 
squared  wave  height.  Group  velocity,  Cg  is  described  by  linear  wave  theory 
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where  C  is  phase  velocity,  d  is  depth,  and  k  is  the  wavenumber  associated  with  peak  fluency  /  of  the 
spectnun. 

Surface  roller  energy  is  given  by  [Svendsen.  1984] 


py4C^ 

21 


(9) 


where  A  is  area  of  the  roller  and  L  is  the  wavelength  of  the  wave.  Engelund  [1981]  uses  an  analogy  with 
a  hydraulic  jump  to  analytically  determine  roller  area, 


A  = 


4dkano 


(10) 


where  a  is  the  angle  of  the  roller/wave  stress  vector,  t,  with  respect  to  horizontal  (Fig  1) .  The  value  of 
a  is  a  function  of  breaker  type  and  is  believed  to  vary  from  a  maximum  value  when  breaking  is  initiated 
to  a  minimum  value  (estimated  to  be  10  degrees)  in  the  inner  surf  zone  [Engelund,  1981].  Thus  Equation 
(9)  becomes 


(II) 

''  8dtano 

The  roller  dissipation  term,  represents  work  done  on  the  roller  per  unit  area  per  unit  time,  and  is 
given  by 
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€,  “  —  *  pgfAtaao  = 

L  4a 

The  work  done  by  the  roller  is  accomplished  by  the  shear  stress  between  the  roller's  bottom  boundary  and 
the  wave’s  surface  boundary,  beginning  at  the  toe  (Fig  1).  Turbulence  generated  at  the  toe  moves 
horizon through  the  top  of  the  wave  and  fills  the  crest  area  as  the  roller  slides  down  the  wave  face 
[Bradshaw  and  Ferris,  1965].  Turbulent  kinetic  energy  is  removed  ftom  the  roller  at  the  same  rate. 

Because  the  generation  of  turbulence  and  energy  dissipation  only  apply  to  breaking  waves,  it  is 
necessary  to  identify  which  waves  are  breaking.  The  average  rate  of  energy  dissipation.  <e,>.  (as  well  as 
the  average  roller  energy  flux,  E.)  is  calculated  by  multiplying  the  dissipation  for  a  single  broken  wave 
of  height  H  by  the  probability  of  wave  breaking  at  each  height,  p*(H),  and  integrating  over  all  H 
[Thornton  and  Guza.  1983].  Breaking  waves  are  identified  by  weighting  the  Rayleigh  function  for  all 
waves,  breaking  and  non-breaking 


p,{H)  =  W{H)p{H)  (13) 

where  fVfH)  is  a  weighting  function,  and  the  area  under  the  distribution  is  equal  to  the  percent  of 
breaking  waves.  The  weighting  function  is  based  empirically  on  acquired  data  which  identify  breaking 
waves,  including  both  incipient  breakers  and  rollers,  and  is  given  by 


W{H)  =  ( 


-=)*l  1 

yd 


cxp(  -  [JLf] 
yd 


(14) 


W(H)<  1  unless  conditions  are  totally  saturated  after  which  all  waves  break  and  W(H)=1. 


Integrating  in  Equations  (11)  and  (12)  through  the  breaking  wave  height  distribution  yields 


<Ht>  = 


<  ^  iydf 


1 


(1  -  ( 


H  ^ 


(15) 


Using  Equations  ( 1 5),  ( 1 1 )  and  (7)  in  Equation  (5)  yields  the  ensemble  averaged  energy  flux  balance 


;  Hi 


d  --mi,  J 

dx  32dtaiio 


1 


yd 


1)  = 


3v^p^  , 

I6d  iyd)^^ 
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I  *  (%)= 
yd 


]) 


(16) 


used  to  describe  wave  height  tiansformation  including  a  surface  roller.  For  saturation  conditions  W(H)=1 
(when  all  waves  are  breaking),  the  ensemble  averaged  breaker  height  reduces  to 


<hI>  = 

4 


(17) 


simplifying  the  energy  flux  balance  equation  to 
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^  a  .3VlipC^C06Cjy3  . 
at  32dimo 


\6d 
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(18) 


For  saturation  conditions,  we  must  determine  when  the  probability  of  all  waves  breaking  is  equal  to  one, 
or 


/  p^(H)d(N)  =  /  Mfmm/f)  =  1 
0  0 


Substituting  Equations  (6)  and  (14)  into  Equation  (19)  and  solving  for  the  shallow  water  depth  limiting 
condition,  saturation  conditions  exist  when 


=  1.27  (20) 
Y<f 

and  Equation  (18)  is  used  instead  of  (16).  A  simple  forward  stepping  numerical  scheme  is  sufficiently 
accurate  [Thornton  and  Guza,  1983]  to  solve  Equations  (16)  and  (18),  written 

«  m.C^,  *  (TO.  (21) 

Due  to  the  equation's  nonlinear  nature,  (H„„,)2  (shoreward)  must  be  solved  iteratively  given  (//„„),,  mean 
incident  angle  a,  and  mean  frequency  f. 
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III.  HELD  DATA 


The  surface  roller  model  is  verified  by  comparing  data  acquired  during  three  field  experiments 
conducted  on  two  near-planar  beaches  of  different  slope  (NSTS)  and  a  barred  beach  (DELILAH).  The 
experiments  used  extensive  alongshore  and  cross-shore  arrays  of  current  meters  and  wave  sensors.  Only 
the  cross-shore  arrays  are  considered  here.  Days  for  data  comparison  are  chosen  from  each  experiment 
based  on  model  assumptions  of  straight  and  parallel  contours,  steady  state  conditions,  and  narrow- 
bandedness  in  frequency  and  direction.  Table  1  summarizes  offshore  wave  and  beach  conditions  for  the 
days  analyzed. 

A  TORREY  PINES  BEACH,  CALIFORNIA 

The  first  NSTS  experiment  was  held  in  November,  19  78  at  Torrey  Pines  Beach,  located  north  of  San 
Diego,  California.  Beach  slope  was  approximately  linear,  varying  between  1:50  in  the  surf  zone  to  1:20 
in  the  swash  region.  Beach  foce  topographies  were  generally  concave-up  and  beach  topography  was 
relatively  uniform  in  the  alongshore  direction. 

A  cross-shore  array  of  1 1  current  meters,  4  pressure  transducers,  and  4  wave  staffs  were  used  to 
measure  rms  wave  heights  across  the  surf  zone.  Dtda  were  telemetered  from  the  instruments  at  Torrey 
Pines  to  the  Shore  Processes  Laboratory,  Scripps  Institution  of  Oceanography  and  digitally  recorded  at 
a  rate  of  64  samples  per  second.  The  data  was  then  low-i)ass  filtered  and  reduced  to  2  samples  per  second 
[Thornton  and  Guza,  1983]. 

A  wide  variety  of  wave  and  weather  conditions  were  encountered  during  the  experiment.  Significant 
offshore  wave  heights,  H,,,  varied  between  .6  and  1.6  meters.  The  average  peak  frequency  of  the  incident 
wave  spectra,  was  feirly  constant  at  0.07  Hz  (Tp=14s).  Incident  wave  angles  in  10-m  depth  were  limited 
to  less  than  IS  degrees  due  to  offrhore  island  shadowing  and  refraction  [Pawka  et  al,  1976].  Two  days, 
4  November  and  10  November,  are  chosen  for  data  comparison. 


9 


B.  SANTA  BARBARA,  CALIFORNIA 

The  second  NSTS  experiment  was  conducted  at  Leadbetter  Beach  in  Santa  Baibaia,  California  over 
a  one  month  period  in  the  winter  of  1980.  Bottom  contours  inside  6m  d^h  were  relatively  straight  and 
parallel.  Mean  nearshore  beach  slope  varied  between  1:33  and  1:16,  depending  on  tide  level. 

A  cross-shore  array  of  16  current  meters  and  6  pressure  sensors  measured  the  cross-shore  wave 
transformation  properties.  Two  wave  slope  arrays  in  deep  water  measured  incident  wave  spectra.  All 
instrument  signals  were  cabled  to  shore  and  digitally  recorded  [Thornton  and  Guza,  1986], 

Leadbetter  Beach  is  oriented  approximately  east-west  on  a  predominately  north-south  California  coast. 
To  reach  Santa  Barbara,  the  open-ocean  North  Pacific  swell  must  pass  between  Point  Conception  and  the 
Channel  Islands,  a  narrow  window  of  ±9  degrees  centered  on  249  degrees  true  [Thornton  and  Guza, 
1986].  This  unique  beach  orierttarion  provided  good  data  sets  3-6  February,  of  which  days  4  and  5  are 
chosen  for  model  comparison. 

C.  DELILAH 

The  third  data  set  used  was  acquired  on  a  barred  beach  during  DELILAH  held  October  1-21,  1990 
at  the  US  Army  Corps  of  Engineer's  Field  Resemoh  Facility  at  Duck,  North  Carolina.  Both  bar  location 
and  beach  slope,  which  varied  between  1:25  and  1:19,  responded  to  a  wide  variety  of  wave  conditions 
[Church  and  Thornton,  1993]. 

Directional  wave  spectra  were  obtained  in  8ra  depth  fiom  an  alongshore  array  of  bottom  mounted 
pressure  sensors.  A  cross-shore  array  of  9  pressure  sensors  were  located  from  4m  depth  to  the  shoreline. 
The  array  provided  near-continuous  data  throughout  the  experiment  with  a  sampling  rate  of  8Hz.  An 
autonomous  Coastal  Research  Amphibious  Buggy  (CRAB)  was  used  for  daily  bathymetric  measurements. 

Three  days  (Oct  10, 12,  and  14)  are  cnosen  for  model  comparison  when  waves  were  relatively  narrow- 
banded,  varying  in  rms  wave  height  from  .77  to  1.30  meters,  and  arriving  at  angles  between  1  and  12 
degrees  from  normal  to  the  beach. 
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IV.  MODEL  RESULTS 

Model  results  for  representative  days  during  each  experiment  are  presented  in  Figures  2*8.  The  two 
free  parameters,  y  and  a,  of  the  roller  model  are  determined  by  fitting  model  to  observed  data  in  a 
least  square  sense.  Optimal  values  for  each  day  and  the  resultant  percent  rms  error  are  given  in  Table  1 . 
Root  mean  square  error  is  determined  at  the  wave  height  observation  locations.  The  predicted  and 
observed  wave  heights,  and  dissipation  as  a  function  of  offshore  distance  and  depth  are  shown  in  Figures 
2a>8a.  The  gradient  wave  and  gradient  roller  energy  fluxes  as  a  function  of  offshore  distance  and  depth 
are  shown  in  Figures  2b-8b.  Dissipation  (modeled  by  Equation  (12))  in  Figures  2a-8a  is  the  sum  of  the 
gradient  energy  fluxes  in  Figures  2b>8b.  The  percent  of  breaking  waves  as  a  function  of  offshore  distance 
and  depth  are  shown  in  Figures  2c-8c. 

The  DELILAH  results  indicate  two  significant  wave-breaking  regions.  The  first  region  just  seaward 
of  the  offshore  bar  is  relatively  broad,  while  the  second  region  shoreward  is  more  abrupt.  Weather  and 
wave  conditions  discussed  earlier  are  evident.  On  October  10,  relatively  small  waves  with  an  average 
period  of  9.1  seconds  arrived  at  an  average  incident  angle  of  1.9  degrees.  The  tide  level  is  .7m.  As  a 
result,  only  55%  of  the  waves  break  at  the  offiihore  bar.  Some  waves  continue  breaking  as  they  propagate 
shoreward,  imtil  at  approximately  20  meters  of&hore  the  remaining  waves  break  and  saturation  occurs. 

Waves  began  building  on  the  1 1th  due  to  a  remote  storm  and  continue  building  through  the  13th.  By 
12  October,  larger  waves  are  approaching  with  an  average  period  of  15.6  seconds  at  incident  angles  up 
to  12  degrees.  The  tide  level  at  observation  time  is  .8  meters.  The  majority  of  waves  (80%)  are  predicted 
breaking  at  the  of^hore  bar  due  to  the  wave  height  increase,  resulting  in  additional  dissipation.  Some 
waves  continue  to  break  shoreward  of  the  bar  until  reaching  saturation  at  approximately  25  meters  offshore 
where  all  waves  break.  Note  how  the  bathymetric  trough  has  widened  by  {q)proximately  20  meters  since 
10  October,  and  beach  slope  has  decreased.  By  October  14,  smaller  waves  are  arriving  at  an  average 
period  of  10  seconds  and  average  incident  angle  of  7  degrees.  At  a  low  tide  level  of  -.2  meters  at 
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observation  time,  saturation  occurs  seaward  of  the  bar.  The  dissipation  region  has  narrowed  and  fewer 
waves  reform  in  the  trough,  uow  widened  by  the  storm  to  SO  meters.  Beach  slope  has  further  decreased 
shoreward  of  the  bar. 

Using  data  acquired  at  5  beaches,  Sallenger  and  Holman,  [1984]  rqrplied  linear  regression  to  conclude 
thm  for  satur^ed  breaking  waves,  y  is  related  to  beach  slope,  p,  (Fig  9)  by  the  relation 

=  3Aanp  +  0.30  (22) 

and  that  ydoes  not  appear  to  depend  on  wave  steepness.  Optimal  y  and  a  values  reflect  wave  and  beach 
conditions  (Table  1)  and  confirm  the  linear  relationship  in  Equation  (22).  How  these  free  parameters 
reflect  changing  conditions  can  be  best  observed  during  DELILAH.  On  October  10th,  beach  slope  is 
steepest  requiring  a  relatively  large  y.  On  12  October,  larger  waves  have  decreased  the  beach  slope, 
reflected  in  the  y  decrease.  Finally,  by  14  October  the  beach  slope  has  decreased  significantly  resulting 
in  a  major  decrease  in  y. 

Roller  steepness,  sincr,  also  changes  with  wave  conditions;  the  steepest  angles  are  required  for 
plunging  breakers,  such  as  those  occurring  on  the  12th,  and  smaller  angles  for  spilling  breakers.  The  same 
trends  hold  for  the  near-planar  Leadbetter  and  Torrey  Pines  beaches,  but  with  smaller  percent  rms  error 
due  to  constant  beach  slopes. 
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V.  DISCUSSION  AND  CONCLUSIONS 


Sensitivity  tests  varying  o  and  y  versus  rms  wave  height  prediction  error  are  shown  in  Figures  10a 
and  lOb,  All  model  runs,  regardless  of  wave  or  beach  conditions,  produce  the  same  trends.  Roller 
steepness  can  vary  up  to  40  degrees  in  most  cases  negligibly  changing  rms  error.  The  reason  for 
insensitivity  to  o  is  seen  by  examining  the  energy  flux  balance  Equation  (16). 

Roller  steepness  is  used  to  describe  both  the  surface  roller  production  and  dissipation  terms.  After 
substituting  tano  into  Equation  (10)  into  the  dissipation  term  Equation  (12),  tana  caiKels  and  remains 
only  in  the  roller  energy  term.  As  shown  in  Figures  2b-8b,  the  gradient  wave  energy  flux  is  more  than 
two  orders  of  magnitude  larger  than  the  gradient  roller  energy  flux  term.  This  effectively  reduces  the 
surftice  roller  model  free  parameters  flam  two  to  one  in  y. 

The  model  is  sensitive  to  changes  in  gamma.  The  expected  range  of  y  using  Equation  (22)  is  .3  to 
.6  for  the  range  of  beach  slopes  from  0-.  1 .  Examples  of  varying  gamma  for  planar  and  barred  beaches  are 
shown  in  Figure  1  lb.  Model  error  is  ^rly  insensitive  for  the  barred  beach  (±2%)  when  gamma  is  within 
the  range  .3  to  .S.  Error  rapidly  increases  for  gamma  outside  this  range.  The  reason  for  the  rate  of  change 
difference  can  be  seen  aiuilyzing  small  and  large  gammas  in  Equation  (16).  As  rms  wave  height  is  raised 
to  the  Sth  power  in  the  roller  energy  and  dissipation  terms  compared  with  the  2nd  power  as  in  the  wave 
energy  term,  a  small  decrease  in  y  causes  the  magnitude  of  the  gradient  roller  energy  flux  and  dissipation 
terms  to  rapidly  approach  wave  energy  flux  magnitude.  This  rapid  change  in  magnitude  causes  rapid 
changes  in  rms  wave  height  prediction-  .^nd  percent  rms  error.  A  large  positive  change  in  y  only  decreases 
the  roller  energy  and  dissipation  terms  more,  which  are  already  small  relative  to  the  wave  term,  resulting 
in  slower  percent  rms  error  rate  changes. 

Relatively  low  values  of  y  compared  to  the  expected  range  (.3-.6)  are  required  to  minimize  the  rms 
error  for  Torrey  Pines  data  (Table  1).  The  reason  is  hypothesized  to  be  linear  bore  theory  deficiencies 
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based  on  conclusions  by  Cacina,  [1989],  A  review  of  the  bore  dissipation  model  is  required  for  further 
investigation. 

The  bore  dissipation  model  [Thornton  and  Ciuza,  1983]  predicts  shoreward  rms  wave  heights  within 
an  rms  error  of  ±9%  for  random  waves  and  arbitrary  bottom  profiles  using  the  same  data  for  Toney  Pines 
and  Leadbetter  Beaches.  The  bore  model  energy  flux  balance  equation  is  given  by  Equation  (3),  D  being 
described  by  linear  bore  dissipation.  Frictional  dissipation  at  the  bottom  boundary  layer  is  shown  to  be 
negligible  (less  than  3%  of  breaking  wave  dissipation  within  the  surf  zone)  except  in  the  run-up  region 
and  it  is  subsequently  neglected. 

Applying  conservation  of  mass  and  momentum  at  regions  of  flow  upstream  and  downstream  of  the 
bore,  the  average  rate  of  energy  dissipation  per  unit  area  is  calculated  [Stoker.  1957] 


(*2  -  h,r 

hA 


Q 


mtQ 


(23) 


where  is  a  breaker  coefficient  of  0(1),  representing  the  fraction  of  foam  on  the  wave  &ce  (accounts  for 
various  breaker  types).  Q  is  the  volume  discharge  per  unit  area  across  the  bore  (Fig  1 1),  described  simplest 
by  Hwang  and  Divoky,  [1970]  as 


(24) 


where  C  is  wave  speed  and  L  wavelength.  Substitution  into  Equation  (23)  yields 
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(25) 


^  pgjBHfC  , 

*  ALh  Ah 

The  bore  dissipation  function  has  the  same  form  as  roller  dissipation  with  the  exception  of  the  B 
coefficient.  The  reason  for  their  similarity  is  because  a  hydraulic  jump  (or  bore)  is  used  to  calculate  the 
area  of  a  roller  (Equation  10)  when  deriving  roller  dissipation.  The  fact  that  they  both  are  similar  after 
derivation  from  two  perspectives  would  lead  one  to  expect  similar  results. 

The  sur&ce  roller  model  differs  from  the  bore  dissipation  model  by  only  the  roller  energy  term. 
Although  the  roller  term  is  much  smaller  than  the  wave  energy  term,  it  is  a  significant  addition  to  the  bore 
model.  The  term  provides  two  changes.  First,  it  allows  the  propagation  of  dissipation  shoreward.  Bore 
model  dissipation  results  for  12  October  (DELILAH)  and  4  February  at  Santa  Barbara  (Fig  12)  are 
compared  with  the  roller  model  (Figs  3a  and  Sa).  The  dissipation  plots  are  similar  for  both  models  at  the 
planar  beach,  and  in  the  offshore  bar  region  at  Duck.  The  significant  difference  occurs  in  the  trough  at 
the  barred  beach.  In  the  bore  model,  waves  stop  breaking  in  the  trough  moving  dissipation  towards  zero. 
In  the  roller  model,  waves  break  at  the  ofEshore  bar  and  the  roller  continues  through  the  trough  and 
nearshore  region,  which  is  mote  in  line  with  observations.  The  net  result  is  the  roller  model  more 
accurately  predicts  shoreward  wave  heights  while  providing  more  accurate  offihore  dissipation  values,  as 
shown  by  the  numerical  scheme  in  Equation  (21). 

Secondly,  the  roller  model  decreases  the  number  of  sensitive  input  parameters.  Thornton  and  Guza, 
[1983]  conclude  simple  bore  theory  underestimates  dissipation,  for  which  the  B  parameter  compensates. 
They  also  conclude  that  B  and  y  could  be  combined  into  one  coefficient,  but  are  left  separate  for  greater 
physical  insight.  Cacina,  [1989]  combined  the  B  and  y  t^nn,  called  the  BG  parameter,  described  by 


15 


(26) 


The  bore  dissipation  model  was  iterated  over  the  BG  parameter  to  reduce  percent  rms  error  in  a  least 
squares  sense.  Results  of  one  barred  beach  day  are  shown  in  Figure  13,  revealing  a  similar  free 
parameter/rms  error  curve  as  the  roller  model. 

Cacina,  [1989]  concluded  that  B,  a  measure  of  wave  breaking  intensity,  is  correlated  with  the  deep 
water  surf  similarity  parameter  described  by  Equation  (1).  It  should  follow  that  the  roller  model's  single 
sensitive  free  parameter,  y,  be  a  function  of  beach  slope  and  the  deep  water  surf  similarity  parameter. 
Results  of  varying  y  with  these  parameters  for  all  days  over  various  bottom  profiles  are  shown  in  Figures 
14-16  and  Table  1.  An  ai^roximate  linear  relationship  exists  between  y,  tan^,  and  C.  confirming  previous 
conclusions  by  Battjes  [1974],  and  Sallenger  and  Holman  [1984].  A  high  correlation  coefficient  of  0.93 
is  found  between  y  and  tan^  (Fig  14).  Gamma  is  not  as  highly  correlated  (.67)  with  wave  steepness 
(fyLo)  (Fig  15)  or  the  deep  water  surf  similarity  parameter  (.58)  (Fig  16) . 

Roller  model  accuracy  may  be  improved  by  using  variable  coefficients  across  the  surf  zone.  Constarrt 
Ys  and  o's  were  used  for  each  run.  Optimal  values  were  chosen  by  fitting  model  generated  H,.,  values 
to  observed  wave  heights.  Improvement  would  require  y  and  a  be  specified  throughout  the  wave 
prediction  region. 

In  applying  the  model  based  on  this  analysis,  using  'f^.31  gives  wave  height  prediction  accuracy 
within  10%  rms  error  for  all  beach  profiles.  For  barred  beaches,  a  value  of  y=.34  is  suggested;  for  planar 
beaches,  use  y=.30  (averages  for  barred  and  planar  beach  days).  a=20  degrees  is  sufficiently  accurate  for 
ail  profiles.  The  lower  than  expected  values  of  y  according  to  Equation  (22)  indicates  that  the  roller  model 
y  is  accounting  for  other  frctors  associated  with  complex  wave  breaking  processes. 
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The  roller  model  improves  the  bore  model  by  decreasing  av^age  rms  error  from  8.6%  to  4.S%  for 
all  breaker  types  and  variable  bathymetry,  accomplished  by  more  accurately  describing  dissipation 
throughout  the  surf  zone.  The  rouer  model  simultaneously  reduces  the  requirement  of  two  free  parameters 
to  essentially  one  in  y. 
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TABLE  I.  ROLLER  MODEL  PARAMETER  RESULTS.  Free  parameter  values  are  given  with 
resulting  predicted  vs  observed  wave  height  rms  error.  1 .  Beach  slope-  mean  slopes  for  plaiuu  beach  days; 
foreshore  slopes  for  barred  beach  days.  2.  Mean  incident  wave  frequency  in  Hertz.  3.  Deep  water  wave 
steepness  parBmeter.4.  Deep  water  surf  similarity  parameter.  S.  Mean  incident  wave  angle  in  degrees.  6. 
The  number  of  wave  height  observations  for  a  given  day.  7.  Roller  steepness  angle-  an  adjustable  free 
parameter  in  the  roller  model.  8.  Gamma-  an  adjustable  free  parameter  in  the  roller  model.  9.  Percent  root 
mean  squared  error  for  observed  vs  model  predicted  wave  heights. 
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C  =  phase  speed  of  wave 
A  =  area  of  roller 
L  =  length  of  jump 
o  =  roller  steepness  angle 
h  =  mean  water  depth 
H  =  height  of  wave/bore 


Figure  I.  Sur^Ke  Roller  on  a  wave  ^e  depicting  stress  vector  (x)  and  roller  steepness  angle  (O) 
with  respect  to  horizontal. 
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Figure  2a.  Roller  model  wave  height  prediction  using  'p.37,  o»30  degrees. 
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Figure  2b.  Roller  model  gradient  wave  and  gradient  roller  energy  flux.  y=.37, 
cr=30  degrees. 
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Figure  2c.  Percent  of  waves  breaking  as  a  ftinction  of  offshore  distance  and  depth. 
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Figure  3«.  Roller  model  wave  height  prediction  using  y=.i4,  o»40  degrees. 
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Figure  3b.  Roller  model  gradient  wave  and  gradient  roller  energy  flux.  *^.34, 
0=40  degrees. 
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Figure  3c.  Percent  of  waves  breaking  as  a  function  of  of&bore  distance  and  depth. 
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Figure  4a.  Roller  model  wave  height  prediction  using  1^.30,  o«30  degrees. 
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Figure  4b.  Roller  model  gradient  wave  and  gradient  roller  energy  flux. 
<J=30  degrees. 


Figure  4c.  Percent  of  waves  breaking  as  a  function  of  offshore  distance  and  depth. 
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Figure  5b.  Roller  model  gradient  wave  and  gradient  roller  energy  flux.  ^^'32, 
0*5  degrees. 
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Figure  Sc.  Percent  of  waves  breaking  as  a  function  of  of&bore  distance  and  depth. 
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Figure  6a.  Roller  model  wave  heigbt  prediction  using  ^*.32,  o«20  degrees. 


Figure  6b.  Roller  model  gradient  wave  and  gradient  roller  energy  flux,  y^.32, 
o»20  degrees. 
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Figure  6c.  Peicent  of  waves  breaking  as  a  function  of  offebore  distance  and  depth. 
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Figure  7a.  Roller  model  wave  bei^  predictkm  using  (r^lO  degrees. 


Figure  7b.  Roller  model  gradient  wave  and  gradient  toller  energy  flux.  ^=28 
0=10  degrees. 


Figure  7c  Percent  of  waves  breaking  as  a  function  of  of&bore  distance  and  dqrth. 
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Figure  8a.  Roller  model  wave  height  prediction  using  ^*.27,  o-lS  degrees. 


Figure  8b.  Roller  model  gradiem  wave  and  gradient  roller  energy  flux.  ']r-27, 
0=15  degrees. 
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Figure  8c.  Percent  of  waves  breaking  as  a  function  of  ofbbore  distance  and  depth. 
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gure  9.  Gamma  and  beach  slope  relationship  {Sallenger  and  Holman 
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Figure  10a.  Gamma  sensitivity  test  using  optimum  0's  from  Table  1. 
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Figure  12.  DELILAH  experiment.  Bore  Model  rms  wave  height  prediction  in  meters  (solid  line  above  zero 
reference  line),  turbulent  kinetic  energy  calculation  or  dissipation  (dashed  line),  zero  reference  (dotted 
line),  and  depth  in  meters  (solid  line  below  zero  reference)  on  days  indicated  (Church  and  Thornton, 
1993). 
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Figure  13.  Variation  in  nns  error  due  to  change  in  BG  parameter  for  a  bore  dissipation  model  nin  on  a 
barred  beach  during  SUPERDUCK  experiment  in  Duck.  N.  Carolina  (from  Cacina,  [1989]). 
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Figure  14.  Conelation  of  y  and  beach  slope,  tanp. 


Figure  1 5.  Correlation  of  y  and  deep  water  wave  steepness  (Hg/L,). 
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Figure  16.  Correlation  of  7  and  the  deep  water  surf  similarity  parameter. 
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